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ABSTRACT * 


A transonic equivalent strip (TES) method ha3 
been further developed for unsteady flow computa- 
tions of arbitrary wing planforms. The TES method 
consists of two consecutive correction steps to a 
given nonlinear code such a3 LTRAN2 ; namely, the 
chordwise mean-flow correction and the spanwise 
phase correction. The computation procedure 
requires direct pressure input from other computed 
or measured data. Otherwise, it does not require 
airfoil shape or grid-generation for given 
planforms. To validate the computed results, four 
swept, tapered wings of various aspect ratios, 
including those with control surfaces, are selected 
as computational examples. Overall trends in 
unsteady pressures are established with those 
obtained by XTRAN3S codes, Isogai’s full potential 
code and measured data by NLR and RAE. In 
comparison with these methods, the TES has achieved 
considerable saving in computer time and reasonable 
accuracy which suggests immediate industrial 
appl ications . 

INTRODUCTION 

Considerable attention has been directed in 
recent years towards the technology development of 
transonic aeroelastic applications with particular 
emphasis on transonic flutter predictions. Based on 
transonic smal 1- d i 3 t ur ba n ce equations (TSDE ) , 
computational methods for unsteady transonic flow 
have been developed extensively both in two and 
three dimensions, notably LTRAN2 (now ATRAN2) (Refs. 
1,2) and XTRAN3S computer codes. Various versions 
of XTRAN3S codes (Refs. 3,4,5) and a number of full 
potential methods (Refs. 6,7) are all in good 
progress. Recently, Guruswamy and Goorjian (Ref. 
3), and Borland and Sotomayer (Refs. 4,8) have 
applied their XTRAN3S codes to Northrop F-5 wing 
planform. Bennett et al . (Ref. 9) have applied 
Langley XTRAN3S code for a RAE swept wing; Ruo and 
Malone (Ref. 10) have applied the 3ame code for LANN 
wing computations. Isogai and Suetsugu (Ref. 6) 
have applied their full potential code to various 
planforms including the AGARD standard RAE wing with 
an oscillating flap. Most of these results have 
shown good agreement with the measured data of NLR 
and RAE. 

Although these methods could produce reasonably 
accurate results, computation efficiency and the 
grid generation procedures remain to be improved 
before these codes can be adopted for industrial 


applications. For flutter predictions and 
aeroelastic optimizations, a more efficient computer 
code capable of rapid computations is sought by the 
aerospace industries, since cost-effectiveness i3 
one of their main concerns. 

Motivated by these considerations, we have set 
forth to develop a simple and more efficient method 
for unsteady, three-dimensional flow computations. 
Our objective is to achieve: (1) computation 

efficiency, (2) flexibility and ease of applications 
for flutter and (3) a unified subsoni c/ transoni c 
method for arbitrary planforms. Consequently, a 
preliminary version of the transonic equivalent 
strip (TES) method has been developed (Ref. 11). 
The results obtained for rectangul ar- wing studies 
show promise for the TES formulation. Hence, it 
prompts further development. In this paper, the TES 
method is further developed so that all procedures 
are automated for aerodynamic computations. Thus it 
can be readily adopted by the flutter prediction and 
the aeroelastic optimization programs. To validate 
the present method, a fairly comprehensive 
comparison with available data is given for a 
number of wing planforms including those with 
control surfaces. 

TRANSONIC EQUIVALENT STRIP (TES) METHOD 

The use of strip concept for unsteady transonic 
computations was first proposed by the ONERA group 
(Ref. 12). A similar strip approach, but involving 
quasi-steady approximations, was recently 
implemented at NLR (Ref. 13). In both cases, 
however, only wing planforms of large aspect-ratio 
are treated and possible applications of their 
methods to low aspect-ratio wings are not forth- 
coming. By contrast, the present TES method has 
developed a more general scheme which could handle 
arbitrary planforms including oscillating control 
surfaces . 

Correction Procedures 

Specifically, the present method consists of the 
applications of two correction steps to a given two 
dimensional code; it could be a nonlinear code such 
as LTRAN2 or it could be a time linearized one. The 
basic correction 3teps are: (1) the mean-flow 

correction applied in the chordwise direction and 
(2) the phase correction in the spanwise direction 
(see Fig. 1 for flow chart). The first correction 
is fully automated by an inverse design procedure 
(IAF2 code) in that the local shock structure is 
properly recovered according to the given mean flow 
input provided by a selected computational method or 
by measured data. In this inverse problem, as 
solved by Fung and Chung (Ref. 14), the velocity 
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potential from integrating the pressure on the slit 
representing an airfoil is known up to an arbitrary 
constant. To determine this constant, a closure 
condition is Imposed, e.g., the resulting slope 
distribution being equivalent to a closed body. 
This constant is being updated during the numerical 
iteration process until a converged solution 
satisfying the closure requirement is obtained. 
Once the slope distribution of the new equivalent 
airfoil is found and the steady flow field fixed, 
unsteady responses can then be calculated by varying 
the slopes to account for unsteady motions. The 
latter step can be accomplished by applying LTRAN2 
code to the equivalent airfoil. 

For the second correction, we make use of the 
three-dimensional linear wave analogy in the sense 
that the phase angles are redistributed along the 
span according to the physical model of acoustic 
wave propagations in a uniform medium. This is to 
say that while the first correction accounts for 
reproducing the nonlinear structure of * the three- 
dimensional mean flow, the second correction is 
responsible for the adjustment of the spanwise phase 
lag of the pressure according to an equivalent 
linear three-dimensional flow. In practice, a 
typical lifting surface method such as Doublet 
Lattice code (Ref. 15) is adopted for the second 
correction. Clearly, shock waves cannot be created 
or destroyed by any process of these corrections. 

We should note that the present terminology of 
’’strip method" is defined only by the stripwise 
computation procedure and is otherwise irrelevant to 
the classical strip theory. The above correction 
procedures clearly indicate that the present TES 
approach is equivalently three-dimensional, since 
there is no restriction to the wing aspect-ratio. 

Justification 

It has been pointed out by Fung and Lambourne 
(Ref. 16,17), among others, that an accurate steady 
state with correct shock jump and location is 
essential for correct unsteady aerodynamic computa- 
tions. It is believed that TSDE in general should be 
adequate for computation of unsteady disturbances, 
which are acoustic signals assumed to be small in 
conventional flutter analysis. However, when 
applying TSDE methods, inaccuracy may occur as a 
result of the local failure at the wing leading edge 
and the limitation in the prediction of the shock 
strength. Since an accurate steady state pressure 
field is desired, an alternative is to find the 
airfoil slopes, or equivalently the airfoil, that 
corresponds to a given pressure distribution. This, 
in turn, suggests the inverse design procedure used 
in the first correction. 

Meanwhile, the unsteady aerodynamics of wing 
flutter at transonic speed is complicated by the 
embedded supersonic region. Disturbances downstream 
cannot be felt directly at an upstream point. As a 
result, the phase lag between unsteady motion of the 
wing and correspond! ng aerodynamic response 
increases as the supersonic region gets larger- 
While this effect is important in the chordwise 
direction of a wing, the existence of supersonic 
region should not affect the way an acoustic signal 
propagates in the spanwise direction. With the 
exception of highly-swept wings, the flow- induct ion 
effect in the spanwise direction should be subsonic 
in nature. It is therefore conceivable that the 
correlation between the aerodynamic responses at 
different spanwise stations is a linear one as 
assumed in conventional flutter analysis, and that 
the unsteady aerodynamic responses due to structural 
deformations of a wing at transonic speeds are 


similar in nature to those at the correspond! ng 
subsonic speeds, except in the atreamwise direction. 

ANALYSIS 


Governing Equations 


The simplest form of the t ime- dependent three- 
dimensional TSDE can be expressed a3 
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where 


A - Mj, k l /6 2 ' 2 , B - M 2 k/6 2 ' 1 , C - K - r$ , 
K ■ ( 1-M^)/5 2 ' 2 , r - (Y+1)M 2 and 
0 - c 2 /b 2 6 2 / *. 


The nondimen3ional quantities and coordinates are 
defined as 


♦ - ♦/(ci 2 '*UJ , 

(x,y,z) - (x/c, y/b, z«* /2 /o), 
t - too and k - uc/U^ , 

where all barred symbols denote the true physical 
quantities, parameters c, b, 6 and o» represent the 
root chord, the semi-span, the airfoil thickness 
ratio and the circular frequency of oscillation, 
respectively. 

The potential p can be split into two 
components, i.e., 


♦ - $ 0 (x,z,t) + <f>! (x,y ,2 ,t; Aa) (2) 


where $ 0 satisfies the nonlinear, two-dimensional 
equation (set D » 0 in Eq. (1)) as can be solved by 
LTRAN2 code; p , is the correction potential 
accounting for the three-dimensional effect 
attributed to a small unsteady disturbance due to 
the amplitude Aa. The t hr ee- dimensional unsteady 
disturbance is assumed to be small as compared to 
the two-dimensional one at all times; hence, the 
nonlinear term in the -equation can be neglected, 
resulting in a linear equation for p l§ i.e. 
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Foregoing physical argument in the TES method allows 
for further simplification of Eq . (3) by ignoring 
the coupling term (4> 0 4>i ) • In so doing, Eq. (3) 
is reduced to the acousti c x e§uation , where <J>, can be 
simply solved by the conventional subsonic lifting 
surface method (Ref. 15). 

To further justify the 3panwise connection 
procedure, it is helpful to investigate the 
characteristic surfaces due to Eq. (1), 


(x .| t) \ (of !)* 2 ♦ (£ * . (p * £>t 2 


(4) 


Expressed in the dimensional form and for t > 0, Eq. 
(4) can be recast into a general form, see Fig. 2A 

(x - U.t) ♦ Uy) ♦ Uz) - Ua w t) (5) 


where A 2 - 1 , for linearized subsonic flow, 

X 2 - WY+DM^M^ , for transonic flow. 
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M is a local, _small-di sturbance Mach num 
defined as M. * <fr-/a w ; hence, X 2 must be posltivfe'fcr 
zero. Eq. (5) indicates that the traveling source 
of small disturbance is emitted from (U^t.0,0) at 
the instant t, whose wave front is a sphere for 
linearized subsonic flow and an ellipsoid for 
transonic flow (Fig. 2A) . Hence, the expanding wave 
fronts propagate at wave speeds in_the chordwise- 
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and P denoting the pressure at upper and 
lower Surfaces , the superscript j - t, N and "8," 
denotes the linear subsonic values and "N" denotes 
the nonlinear transonic values. Clearly, 
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In terms of the strip concept, Eq , 
written as 


Ac 


Pi 
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c, P*i (x ,y ,z; k ) 
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(6) can be 


(9) 


( 10 ) 


Note that at x * U^t or y. - 0, the rate of 
propagation becomes identically sonic; otherwise it 
is time-dependent. Most importantly, along these 
strip planes, the event diagrams (x-t diagrams) 
clearly show that while the propagation in the 
chordwise strip plane is of mixed type which could 
be either locally subsonic or supersonic, depending 
on the difference between a^ and U^, the propagation 
in the spanwise strip is unconditionally subsonic, 
which is independent of the a and U relations 
(Fig. 2B) . 

Boundary Conditions 

On the mean surface of the wing planform the 
potentials must satisfy the tangency condition, i.e. 
at z - 0, 


where c, and c 0 are constants and the subscripts 0 
and 1 denote, respectively, the two-dimensional 
strip value and the three-dimensional value at the 
spanwise location y - y.. While the first 
correction is applied at tVie level of Eq. (8), the 
spanwise correction is implied by the pressure-mode 
relation between pi and po which can be expressed as- 

Pi * Po f J (y;k) . (11) 

According to the foregoing arguments of wave 
propagations in the stripwi3e characteristic plane, 
the nonlinear spanwise pressure function f (y;^<) can 
be approximated by its linear counterpart f (y;k) 
throughout the spanwise correction procedure. 

In passing, we note that in all the figures 
presented, Ac * and Ac n are the real and imaginary 
parts of the unsteady pressure Ac representing the 
in-phase and out- of- phase pressure coefficients. 


■ 

4F x (x,y i ) + a, ♦ a, •(H )[ + kHj.) , 

(6) 

*‘z ' 

ia • (H x * kH fc ) , 

(7) 


where F(x,y ) is the stripwise wing surface 
geometry, y^is the "ith" spanwise location, a 0 is 
the mean ang^e of attack, and H(x,t) depicts the 
wing motion with oscillation amplitude a, ♦ Aa. 
While Aa is the small amplitude which induces the 
three-dimensional unsteady disturbances, Aa is 
actually related to the two-dimensional amplitude a x 
by Aa * o(a! ) . 

Outside of the wing planform, the potential 4> 
must satisfy the radiation condition in the lateral 
and the upstream far fields. The zero pressure jump 
condition across the wake sheet and along which 
there is no flow discontinuity must be maintained. 
Also, an unsteady pressure wave must attenuate far 
downstream. Since these boundary conditions are 
linearized consistent with the small-disturbance 
assumption, expressions for $ 0 and t become 
decoupled. Thus, it renders the correction steps to 
be applied in a consecutive manner. 

Pressure Coefficients 

Due to a small oscillatory amplitude a, the 
pressure coefficient can be decomposed into 

c - c + A c • a 
P P P 

where c is the steady mean pressure coefficient and 
A^ is Ehe unsteady pressure coefficient defined as 


RESULTS AND DISCUSSION 

To demonstrate the present TES method and to 
validate its computed results, four different wing 
planforms are selected for comparison with available 
data. These include Northrop F-5 Wing in pitching 
motion and with an oscillating flap, LANN wing in 
pitching motion, RAE/AGARD tailplane in pitching 
motion and AGARD standard RAE wing with an 
oscillating flap. The sections selected are given 
in the table 1 below; 


Table 1 Selected Wing Sections 


Planforms 

Sections 

/ 

Semi-Span % 


F-5 Wing 

1, 2, 3, 5, 

/ 

18.1, 

35.2, 

51.2, 

72.1, 


6, 7, 8 

/ 

81.7, 

87.5, 

97.7 


LANN Wing 

2. 3, K. 

/ 

32.5, 

47.5, 

65.0 



5, 6 

/ 

8a. 5 . 

95.0 



RAE Tailplane 

1. 2, 3, 5 

/ 

Ilt.o, 

42.0, 

66.0, 

96.0 


Steady Mean-Flow Results : Equivalent Airfoil Design 

For all steady mean-flow pressure design 
outputs, the present computed results using the IAF2 
code for equivalent airfoil design are shown by 
solid lines, whereas the input data in general are 


142 



ORIGINAL PAGE IS 
OF POOR QUALITY 


shown by open symbols. Figure 3A and 3B display the 
pressure inputs for upper and lower surfaces at 
selected sections listed above at H m ■ 0.9 for LANN 
wing and F-5 wing based on NLR* s "measured data 
(Refs. 18,19). For the case of F-5 wing at M b - 
0.95, it can be seen from Figure 3C that "strong" 
shock occurs near the wing trailing edge. Hence 
some care must be exercised to input the pressure 
data. To demonstrate the flexibility of the TES 
method, we adopt the computed results of XTRAN3S- 
Ames (Ref. 3) as our pressure inputs. Notice that 
minor discrepancies between input and output results 
appear behind the mean shock. Figure 10A shows the 
comparison of input data measured by RAE (Ref. 9) 
and the TES output results for a RAE tail plane at 

- 0.9. Figure 11 presents the steady pressures 
for the RAE wing at 45* spanwise location using the 
computed results of the Bai 1 ey-Bal 1 haus code (or 
GACBOPPE code) as inputs, as no measured data were 
provided for this section in Ref. 20. For all cases 
considered, the comparison between the input data 
and the output computed results are generally in 
good agreement . 

Unsteady Pressure Results? LTRAN2/TES Computations 

For unsteady computations, the LTRAN2 code Is 
adopted as our computation basis because of its 
inclusion of nonlinearity and ease of application. 
In all figures presented for unsteady pressures, 
solid lines denote the present TES method, whereas 
lines attached with triangular symbols represent 
various versions of XTRAN3S codes, with the 
exception of Figure 11. The open square and circle 
symbols denote the NLR or RAE measured data for in- 
phase and out- of- phase pressures, respectively. 

Figures A, 5 and 6 contain unsteady results at 
various spanwise sections of Northrop F-5 wing in 
pitching oscillation at M w • 0.9 and at given 
reduced frequencies k - 0.27*1, 0.55 and 0.136, 
respectively. Figure 7 presents the unsteady results 
for the same wing in pitching oscillation at M^ - 
0.95. The pitching axes for all cases of F 75 wing 
are located at 50* root chord. For all cases at M # 

- 0.9, it can be observed that the present results 
practically follow the same trend as those of 
XTRAN3S codes. Better correlations are found at 
lower frequencies than for those at higher 
frequencies. In Figure 5, o verpredi cted unsteady 
shock appears in section 5 which suggests that the 
three-dimensional, wave cancelling mechanism acts at 
high frequencies more effectively than a locally 
two-dimensional one. In Figure 7, the resulting 
unsteady pressures at M m - 0.95 appear also to 
follow the same trend as those of NLR and XTRAN3S 
code, except that the unsteady shock strength is 
again o ver pr edi ct ed by the present method. To 
investigate this problem firther, we change over the 
pressure input as predicted by XTRAN3S to the NLR 
measured data which contains a weaker mean shock. 
As shown in Figure 8, the unsteady pressures are 
indeed improved near the shock and are otherwise 
unaffected. This further verifies ou* contention 
that the steady shock strength and position are the 
most crucial for unsteady pressure predictions. 

Figure 9 presents the in-phase and out- of- phase 
pressures of the LANN wing with pitching axis 
located at 62* root chord. Throughout five spanwise 
locations considered, the present results for upper- 
surface compare more favorably with the NLR measured 
data than do the XTRAN3S results. Meanwhile, 
subcritical flows are predicted for lower surfaces; 
hence, the unsteady pressures do not contain shock- 
jump. Figure 10 presents pressure results for a 


highly-swept RAE tailplane of A - 50.2° at M w = 0.9 
and at zero mean incidence. Because the present 
method uses the measured data inputs, it can be seen 
that the predicted unsteady shock positions 
correlate better with the RAE data than do the 
XTRAN3S results. It should be cautioned that Eq. 
(1) may not be suitable for wings with large 
sweepback angles, as the side wash could be of the 
same order as the mean convective velocity. One 
would therefore expect that for both steady and 
unsteady flow, discrepancies will occur near the 
wing tip, particularly for wings with low aspect 
ratios. Such det er i or at i on in pressure estimates 
was observed in Ref. 9. 

In general, overall trends of unsteady pressures 
were obtained with TES method from Figures 4 to 10. 
In some cases, insufficient adjustment of the phase 
angle causes the under pr edi ct i on of the unsteady 
pressure level. It is believed that this type of 
discrepancy may result from the linear approximation 
inherent in the spanwise correction procedure. 

Oscillating Flap 

Figure 11 presents the computed and measured 
data at 45* semi-span of AGARD standard RAE wing 
with an oscillating flap starting at 70* chord. It 
is interesting to observe that the present results 
and Isogai's full potential method (Ref. 6) are in 
‘good agreement with RAE measured data. Figure 12 
compares in- and out -of- phase pressures at two flap 
sections of F-5 wing at M # - 0.9; the hinge line is 
located at 82* chord. Closed agreements are found 
with XTRAN3S results of Sotomayer and Borland (Ref. 
8) and NLR measured data (Ref. 21). It should be 
pointed out that in presenting the flap oscillation 
data, usual practice in computational methods is to 
connect the data points across the hinge line, in 
the same manner as connecting the shock points; such 
are the cases of Refs. 6 and 8. However, control 
surface singularities in subsonic flow have been 
well established by White and Landahl (Ref. 22). 
Since the flow is assumed inviscid and is locally 
subsonic, hinge line singularity, as we have 
presented in these figures, should prevail. 

Computation Time 

In performing the first correction of the 
present TES code (LTRAN2/IAF2 ) , we used typically 
103 x 97 grid points and assigned 240 time steps for 
each cycle. Pressure data are read usually in the 
fourth cycle, as the aerodynamic response normally 
becomes periodic or harmonic after the third cycle. 
Typically, it takes 700 iterations to achieve the 
IAF2 steady pressure output. In an IBM 3081, the 
CPU time for the first correction amounts to some 
620 seconds. With 10 x 10 panels, the Doublet 
Lattice code requires 40 seconds, with the 
additional procedure, the CPU time amounts to 
roughly 100 seconds. Hence, the total CPU time 
required for running TES code on one strip adds up 
to 720 seconds. The same case (e.g., Northrop F-5 
wing) computed by XTRAN3S code would require 
reportedly 2000 to 4000 seconds CPU time in the 
CRAY-1S supercomputer. This would amount to an 
equivalent CPU time in an IBM 3081 of about 10 hours 
or more. With four strips chosen in this case, 
present TES method usually takes up no more than 
3000 seconds. Therefore, a saving of about ten- to 
twelve-fold in CPU time can be achieved. With 
further improvement to the TES pilot code, it is 
expected that at least another reduction factor of 
two in CPU time can be achieved. 


143 



CONCLUDING REMARKS 

A TES method has been developed for unsteady 
transonic computations about wing planforms. Four 
swept, tapered wings of various aspect ratios 
including control surfaces are selected as 
computational examples for validation of the TES 
computer code. Computed results of TES code 
practically follow the same trends as those of 
XTRAN3S codes and f ul 1- potenti al code; all of them 
have been verified with available measured data by 
NLR and RAE. 

In view of the satisfactory results obtained and 
the effective procedures established by the TES 
method, we believe that it is useful for immediate 
aeroelastic applications. To summarize, the 
following special features of the TES method are 
worthy of notice: 

Applicability to General Planforms 

In addition to its applicability to rectangular 
wings, the TES method is equally applicable to 
swept and tapered wing planforms with any given 
aspect-ratio, including those with control surfaces. 

No Need for Grid Generation 

Any time-domain three dimensional computational 
methods generally requires a grid generation 
procedure, which could be planform- dependent in most 
cases. The present TES method does not require such 
a procedure. 

Computation Efficiency 

A rough estimate in CPU time indicates that to 
compute aerodynamics for one given mode, using the 
present TES code, is at least ten times faster than 
using the XTRAN3S code. With f urther . i mpr ovement of 
the current TES pilot code, it is expected that at 
least another reduction factor of two in CPU time 
can be achieved. 

Flexibility and Ease of Application 

Unlike other unsteady computational methods, TES 
method makes use of the steady-flow pressures 
supplied either by measurement or by steady 
computational method. The flexibility of the TES 
method lies in the pressure input scheme which does 
not require airfoil shapes. For ease of applica- 
tion, the input format of TES code will be unified 
with that of the subsonic doublet lattice code. 

Transonic A.I.C, and Flutter 

With the exclusion of chordwise bending modes, 
the present TES method can be extended to the 
construction of a three-dimensional aerodynamic- 
influence-coefficient matrix, hence, the generalized 
forces. These are the essential building blocks for 
the TES method to become an efficient aerodynamic 
tool for flutter analysis. 
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(A) LANN WING 
M = 0.82 


(B) NORTHROP F-5 WING 
M = 0.90 


(C) NORTHROP F-5 WING 
M = 0.95 


Steady Pressure Inputs and Equivalent-Airfoil Outputs at Various Spanwise Locations for: 

(A) LANN WING at Mean Incidence a 0 - 0.62° 

(B) and" (C) NORTHROP F-5 WING at Mean Incidence cu - 0° 

(□ Upper and A Lower Surfaces - Measured Data: (A) NLR/Horsten, (B) NLR/Tijdeman and 

(C ) Computed Data XTRAN3S/Ames; Present Equivalent-Airfoil Output) 
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SEC 1 SEC 5 SEC 7 



x/c x/c x/t 


Fig. 5 NORTHROP F-5 WING Comparison of In-phase and Out-of-phase Pressures on Upper Surface at Three 
Spanwise Locations: Pitching Oscillation about Mid Root-Chord at Mach Number * 0.90 and 
Reduced Frequency k = 0.550. 

( Present TE§ Method; — A— XTRAN3S/Ames; □ Real and 0 Imaginary - NLR/Tijdeman et al . 

Measured Data) 



Fig. 6 NORTHROP F-5 WING Comparison of In-phase and Out-of-phase Pressures on Upper Surface at Three 
Spanwise Locations: Pitching Oscillation about Mid Root-Chord at Mach Number M # - 0.90 and 
Reduced Frequency k - 0.136. 

( Present TES Method; Q Real and © Imaginary - NLR/Tijdeman et al . Measured Data) 
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FiR. 7 NORTHROP F-5 WING Comparison of In-phase and Out-of-phase Pressures on Upper Surface at Three 
Spanwise Locations: Pitching Oscillation about Mid Root-Chord at Mach Number M w - 0.95 and 

Reduced Frequency k - 0.528. 

( Present TES Method; A — XTRAN3S/Ames ; 0 Real and © Imaginary - NLR/Tijdeman et al . 

Measured Data) 




(A) XTRAN3S/Ames Input (b) NLR Measured Data Input 




Fig. 8 NORTHROP F-5 WING Comparison of Steady Pressure Distributions, In-phase and Out-of-phase 

Pressures at 1 8 % Spanwise Locations: Pitching Oscillation about Mid Root-Chord at Mach Number 

M - 0.95 and Reduced Frequency k - 0.528. 

( - Present TES Method; A- - XTRAN3S/Ames ; 0 Real and © Imaginary - NLR/Tijdeman et al . 

Measured Data) 
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(B) LANN WING Lower Surface 


FiS* 9 LANN WING Comparison of In-phase and Out^of-phase Pressures at Five Spanwise Locations ; 

Pitching Oscillation about 62 % Root-Chord at Mach Number M * 0.82 and Reduced Frequency 
k Q - 0.205 

( Present TES Method; — A— XTRAN3S/Malone; PH Real and Q Imaginary - NLR/Horsten et al 

Measured Data) 





Fig. 


(B) LANN WING Lower Surface [Continued] 
9 LANN WING [Continued] 
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(A) Steady Pressure Distribution 



(B) Pressure Magnitude 


o 

o 



ID STANDARD RAE WING WITH OSCILLATING FLAP 

Steady Pressure Inputs and Equivalent-Airfoil Pressure Outputs at 45J Spanwise Locatic 
Mach Number = 0.90 and Mean Incidence a 0 - 0° 

( □ GACBOPPE*Code Input; Equivalent-Airfoil Output) 

Pressure Magnitude and 

Phase Angle on Upper Surface of Spanwise Location for AGARD RAE Wing with Oscillat 
Flap at Mach Number M^ * 0.90 and Reduced Frequency k * 0.705 with Flap Amplitude 5 = 
( Present TES Method; Isogai's Full Potential Code; o + RAE Measured 
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x/c x/c 

(A) 18% Spanwise Location (Section 1) 




X/C x/c 


(B) 51% Spanwise Location (Section 3) 


Fig. 12 NORTHROP F-5 WING WITH OSCILLATING FLAP 

Comparison of In-phase and Out-of- phase Pressures with Hinge Line at 82 % Chord at 
Sections 1 and 3. 

( Present TES Method; A — XTRAN3S/Sotomeyer; [T| Real and 0 Imaginary 

- NLR/TiJdeman et al . Measured Data) 
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